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Abstract 

The effects of the highly damped modes in the energy and reaction rates 
in a plasma are discussed. These modes, with wavenumbers k 3> kjj, even 
being only weakly excited, with less than fceT per mode, make a significant 
contribution to the energy and screening in a plasma. When the de Broglie 
wavelength is much less than the distance of closest approach of thermal 
electrons, a classical analysis of the plasma can be made. It is assumed, in the 
classical analysis, with h — > 0, that the energy of the fluctuations Tito <C fefiT. 
Using the fluctuation- dissipation theorem, the spectra of fluctuations with 
ft ^ is appreciably decreased. The decrease is mainly for the highly damped 
modes at high frequencies (~ 0.5 — 3ksT). Reaction rates are enhanced in a 
plasma due to the screening of the reacting ions. This is taken into account by 
the Salpeter factor, which assumes slow motion for the ions. The implication 
of including the highly damped modes (with fi ^ 0) in the nuclear reaction 
rates in a plasma is discussed. Finally, the investigations presently done on 
these effects in particle simulations with the sheet model and the multiparticle 
quantum simulation code are described. 



I. INTRODUCTION 



Plasma effects can drastically modify important phenomena in astronomy such as the 
energy density in the early universe and the neutrino flux produced in the interior of the 
Sun. Plasma physics has been extensively studied in the laboratory. However, in astronomy, 
research on plasmas has received little attention. The lack of a more consistent investigation 
of the plasma effects is particularly true for plasma effects in basic processes such as reaction 
rates and equation of state. As all stars possess very hot and dense plasmas, plasma effects 
on reaction rates and equations of state can significantly affect stellar evolution calculations. 
These effects can also influence the early universe, especially at the radiation epoch, when 
the amount of energy in the plasma is comparable to the amount of energy in the radiation 
(photons). 

We discuss in the present article the effects of quasi modes and fluctuations in plasma on 
energy and reaction rates. Even non-magnetized plasmas contain an oscillating magnetic and 
electric fields that produce fluctuations in density for example. The fluctuation-dissipation 
theorem is a powerful description of the fluctuations in plasma for it naturally includes all 
the fluctuations, such as the quasi-modes. 

What are quasi modes or highly damped modes? This denomination refers to modes 
that are not the normal modes of the plasma. For example, in the case of a non-magnetized 
plasma, all modes that are not plasmons (in the longitudinal component) nor photons (in 
the transverse component). We would expect that the energy in the fluctuations resides 
mainly in the normal modes, around the plasma frequency for example (in the case of 
longitudinal component of a non magnetized plasma) and that the contribution of the quasi- 
modes would be negligible. We show that the contribution for the Debye-Huckel correlation 
energy for instance, is appreciable. This is crucial since by not taking the limit of H — > 
for classical plasmas, the fluctuations carried by the highly damped modes are strongly 
affected. Therefore, by decreasing the energy in the highly damped modes, the energy of 
the plasma is strongly affected. The fluctuations in plasma can affect also the screening and 
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the enhancement of nuclear reaction rates. We will comment on that and on the role of the 
highly damped modes. 

Presently, we are investigating the effects of fluctuations in particle simulations using two 
approaches: through a classical description with the sheet model ]]]] and through a quantum 
description with the multiparticle quantum code 0. The plasmas in the interior of stars are 
dense enough that quantum effects are important. Many analytical calculations to include 
quantum effects have been attempted. However, all these calculations suffer from several 
limitations, such as not including for example, non-linear effects. Although many studies 
have used classical codes to investigate plasma, currently there are no existing quantum 
particle codes that can directly simulate basic plasma effects [[|. The only existing code 
deals with very few particles (~ 2). The multiparticle quantum code is a cutting edge code 
to investigate quantum effects in plasma being able to handle hundreds of quantum particles 
and analyze a wide spectrum of frequencies and wavenumbers. 

The outline of this paper is the following: in the Section II, we describe the electromag- 
netic fluctuations; in Section III and IV, we comment on the effect of the highly damped 
modes on the energy of a plasma and in reaction rates. We will concentrate on the recent 
work of Opher and Opher (1999) || that used a semi-classical treatment; in Section V, we 
describe our present investigation of those effects with the sheet model and the multiparti- 
cle quantum code; and finally, in Section VII, we draw the conclusions commenting on the 
implication of those effects in astronomy and in the early universe. 

II. ELECTROMAGNETIC FLUCTUATIONS 

Plasma presents electromagnetic fluctuations even in thermal equilibrium and in a non- 
magnetized state. Detailed study of the electromagnetic fluctuations has been performed 
by Dawson ||, Rostoker et al. 0, Sitenko et al. [[/J, and Akhiezer et al. ||. Most of these 
results are compiled in the books of Sitenko and Akhiezer et al. [pyT0|. There has been a long 
tradition in plasma studies to use the analysis of the fluctuations in a plasma as a powerful 
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tool to investigate intrinsic properties of the plasma like the amount of energy, screening 



and diffusion. For example, the studies of Ichimaru 11 1, ONeil and Rostoker [12|, Lie and 



Ichikawa ]13| of the correlation energy of a plasma; of Dawson and Lin et al. [II] on the 



energy and screening properties, and the analysis of anomalous electron diffusion using the 
magnetic fluctuations. The fundamental approach to determine the fluctuation properties 
of any system in the linear regime is the fluctuation-dissipation theorem. 

The level of thermal fluctuations in a plasma is related to the dissipative characteristics 
and temperature of the plasma. Plasma waves are emitted by particles as they move about 
the plasma and are absorbed (e.g., Landau damped) by the plasma. The balance between 
emission and absorption leads to a thermal level of field fluctuations. The fluctuation- 
dissipation theorem describes this level of fluctuations that includes, besides the normal 
modes of the plasma, modes that have a short lifetime, transient modes. For short lived 
modes as the damping is faster, the excitation has to be faster so that the emission is equal 
to the absorption in thermal equilibirium. 

The quasi-modes (or heavily damped modes) are modes with high wavenumber k ^> kp. 
Therefore, we are talking about physical phenomena that exist inside the Debye sphere. 
For the plasmas treated, the closest distance between the particles d is much less than the 
Debye sphere, i.e., k max = l/d > ko- For example, for T = 10 5 K, n = 10 13 cm -3 , 
kmax = 1-1 x 10 4 /cd. For these plasmas, the number of particles inside the Debye sphere is 
large ~ 10 4 (g ~ 10~ 4 ). The particles travelling inside the Debye sphere feel the fluctuating 
electric fields of the bare particles and are continuously scattered by these fields. As they 
move, they emit plasma waves. These waves have a very small phase velocity (as k 3> kp), 
less than the thermal velocity of the other particles and are heavily Landau damped. By 
suffering Landau damping they transfer energy to other modes. The particles excite new 
modes which are again short lived and the process is repeated. As noted by Tajima et al. 
[15], "an individual mode (or field) decays by a certain dissipation, giving up energy to 
particles or other modes, while particles (or other modes) excite new modes and repeat the 
process and the amount of fluctuations is related to the dissipation." 



The fluctuation-dissipation theorem gives for the spectra of the electromagnetic fluctu- 
ations in an isotropic plasma in thermal equilibrium 0, 
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where and £<r are, respectively, the longitudinal and transverse permittivities of the 
plasma. The first and second term of Eq. (1) are the longitudinal and transverse electric 
field fluctuations, respectively. The advantage of using the fluctuation-dissipation theorem is 
that we are able to estimate the energy in the electromagnetic fluctuations for all frequencies 
and wavenumbers. This estimation requires the dielectric permittivities obtained from the 
description of the plasma. The calculation includes not only the energy in the fluctuations 
that are normal modes of the plasma, such as plasmons and photons, but also the energy in 
the fluctuations that do not propagate, or are short-lived, i.e., quasi-modes. The expressions 
given in Eqs. (l)-(2) are a quantum description of the plasma 0. Works such as Opher 
and Opher || are semi-classical descriptions because even though the expression for the 
fluctuations was quantic, the dielectric permittivities used were classical. For a relativistic 
or quantum description, it is necessary to obtain adequate dielectric permittivities. This is 
true also if we are interested in including non-linear effects in the dielectric permittivities. 



Recently, Tajima and collaborators revisited the study of fluctuations. They dis- 
cussed a very interesting suggestion that the magnetic fluctuations could result in a spon- 
taneous magnetic field in a plasma in thermal equilibrium since the magnetic spectrum in 
a plasma at low frequencies (~ 0) has an intense peak. The plasma description used was a 
cold plasma description with the collisions described by binary interactions. A later descrip- 
tion by Opher and Opher (1997a) [[HJ including thermal effects (in addition to collisions), 
showed that in fact, the magnetic spectrum, is not a blackbody spectrum being distorted 
at low frequencies. However, the high intensity peak at low frequencies is strongly damped. 



Extending this study to the transverse electric field, Opher and Opher (1997b) |L7] showed 
that the electromagnetic spectrum (the sum of the magnetic and the transverse spectrum) 
is different from a blackbody spectrum in vacuum. It has a plateau at low frequencies, 
possessing more energy than the blackbody spectrum in vacuum. This shows the power of 
a tool like Fluctuation-Dissipation Theorem. 

III. ON THE ENERGY AND SPECTRA OF THE FLUCTUATIONS FOR A 

CLASSICAL PLASMA WITH h + 

The classical theory states that the energy of a plasma in thermal equilibrium is given 
by the ideal gas thermal energy, 3/2 fcgT plus a correction due to the interactions of the 
particles in a plasma, E corr . To examine the question of the energy of a plasma through the 
fluctuations we need to estimate the energy contained in the fluctuations. Therefore, one 
must integrate Eqs. (1) and (2) over wavenumber and frequency, with a specific dielectric 
permittivity. This will give us the energy in the electric field (longitudinal and transverse) 
and in the magnetic field. The quantity that is interesting, in fact, is not the energy in the 
fluctuations, but how much the fluctuations contribute to the plasma energy. For that, it 
is necessary to subtract from (Ef) the energy of the particles when they are far apart from 
each other. This give us what is called the correlation energy, E corr , of the particles. The 
energy of the particles that are separated from each other, at infinity, E in f, is given by 

E mf = -J dk— . (3) 

Subtracting E in f from the energy in the fluctuations, (Ef) give us the correlation energy, 

Ecorr = (Ef) — E in f ■ (4) 



The standard classical description [ID.I2]] uses the limit of % — > implying that huo <C ksT. 



In this limit Eqs. (1) and (2) turn out to be proportional to (ksT/uj) and can be integrated 
in wavenumber in a general way with the aid of the Kramers- Kronig relations without even 
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specifying the particular dielectric permittivity. For instance, the longitudinal component 
turn out to be, 



n2x nk B T 1 1 l 



In the case of a collisionless electron plasma, 

where g is the plasma parameter. For k 3> ko, we obtain the standard result for the 
correlation energy, the Debye-Huckel energy (for higher orders in g see [[l2|j). 

E corr = ~nT (—g) ■ (7) 



2 U27T 

Opher and Opher (1999) ||] recently obtained the interesting result that the energy in 
a plasma is orders of magnitude higher than the Debye-Huckel result. For example, for 
T = 10 5 K and n = 10 13 cm -3 , E corr obtained was 10 5 larger than (p/127r)(3/2nT). They 
calculated in the same way decribed above, the correlation energy for classical plasmas but 
assumed a finite H. That is, they used the general expressions of the fluctuation-dissipation 
which is a quantum description for the plasma but the dielectric permittivity was a classical 
one. Therefore, the description of the particles, electrons and ions was classical, but the 
fields (photons) were quantized. The plasmas studied were classical in the sense that the de 
Broglie wavelength was much smaller than the closest approach between particles in plasma. 
In the purely classical case, when h — > 0, the Debye-Huckel theory is recovered. This study 
showed that, quantum corrections even to classical plasma, can be very important. 

In this section we summarize our recent article on reviewing this work [pl| trying to 



explain in a simple way, why this result was obtained and what was exactly the role of the 
highly damped modes. Some of the questions raised by this result were: if this extra energy 
was due to the cut-off used in the wavenumber integration of the fluctuations; what was the 
effect of a non-zero h; and which modes were responsible for this extra energy: the highly 
damped modes, such as modes higher that the peak of the blackbody, where Hu > ksT. 
Below we summarize our recent analysis 111 



(a) The first point was to clarify the role of the highly damped modes in classical plasmas 
and about the cut-off used in the wave number integration. When integrating the spectra, 
a cut-off is necessary since the fluid description breaks down at small distances. In fact, we 
show that the extra energy is not due to the cut-off in wavenumber. At first, it might seem 
confusing given that as works like Tajima et al. [[13] took it as being the Debye wavenumber. 
The blackbody spectrum is only obtained at high frequencies if we allow the wavenumber 
to be much higher than the Debye wavenumber. In the case of Tajima et al. fig] they 
had to impose a blackbody spectrum at high frequencies and fit it with the spectra at low 
frequencies. Opher and Opher || used k max being one over the closest distance between 
the particles in the plasma and obtained the blackbody naturally. However, by examining 
Eq. (||) it is clear that the Debye-Hiickel theory already assumes that k ^> hp. Therefore, 
even in the pure classical case with fi — > 0, the highly damped modes are included in the 
standard treatment and contribute to the energy. This can be seen in Fig. 1 that shows the 
longitudinal spectra (Ef)^/ (Ef)^^ for the classical case with h — > [Eq. (1)] vs. u)/oj pe for 
diferent k's. The parameters are T = 10 5 K, and n = 10 18 cm -3 and peak of the blackbody 
is 232u pe . The full curve is for k = COl/c^ , the dotted for O.lkp, the dashed for lOkp and 
the long-dash is for k = 100k D . Since as plasmons occupy the phase space for uj ~ u pe and 
for k <C ko, we can see the great contribution of the quasi- modes extending to frequency 
around the peak of the blackbody. The knee in all the curves is due to the Landau damping 
at uj ~ kv te - As expected, the knee moves to the right for higher k's. This emphasizes that 
even for a classical calculation, with 7i — > 0, the modes with frequency different from the 
plasma frequency, carry a substantial portion of the energy. 

(b) What is the effect of not taking the limit of h — > 0? In Eq. (1) we can see that 
the presence of h affects the factor h/ (exp(hu) / kgT) — 1) (since the permittivity used is 
classical). This factor in the limit of h — > turn out to be ksT/u. In Fig. 2 we plotted the 
spectrum (Ef) w / '(El) (t)=0 vs. uj/u pea k, where uj pea k is the frequency peak of the blackbody 
for k = k max , in the cases when with h — > (full curve) and with % ^ (dashed curve). 
We can see that even at frequencies less than the peak of the blackbody this factor has a 
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major influence. The spectra with finite % substantially decreased from the one with % — > 0. 
This shows clearly that quantum effects are important even for classical plasmas at classical 
frequencies {Tiuj <C ksT). 

(c) Why Opher and Opher [|J obtained this extra energy? In Fig. 3 we plotted the 
longitudinal energy, (E^)^ , as a function of the maximum frequency of the integration, 
(already integrated in wavenumber). Both the cases with h — > (full curve) and h ^ 
(dashed curve) are plotted. As we increase the maximum frequency of integration, the 
energy increases as expected, approaching a plateau since we ended up summing all the 
modes in the plasma that carry the energy. The dotted curve is the energy of the particles 
separated at infinity, E in f calculated classically [Eq.@]. The difference between the dotted 
curve and the full and dashed curve is the correlation energy (Eq. (|4])) and it is as expected, 
a small value cx g, the Debye-Huckel result. The energy spectrum in the case of % — > 
approaches the dotted curve, E in f while the case with h ^ does not approach the dotted 
curve but a different plateau value. This is why Opher and Opher [|J obtained an energy 
so different from the standard calculation. The great difference in energy with % — > and 
the case with h ^ is concentrated in frequencies from 0.5 — 2.Qu pea k- Also, this imply 
that the energy of the fluctuations with finite h ^ is in fact reduced. So where is the 
extra energy coming from? The answer is in the energy of the particles at infinity. We 
have to clarify that if by not taking the limit of h —>■ to be finite, the energy of the 
particles at infinity is still given by expression Eq. (|3|). We expect that when the plasma 
parameter g (or the density) approaches zero, A = {^nT)E corr will tend to zero as there 
are no interactions. In fact lowering g, A approaches zero, where A is the Debye-Huckel 
result. However, the same does not occur with A, where instead it approaches a constant 
value ~ 0.06 approximately the plateau approached by the dotted curve. The energy of a 
plasma is given by E = |nT + |A = E^ + ((E 2 ) — E in f), where Eid is the energy of an ideal 
gas, and Ei n f the energy of the particles at infinity. The fact that (E 2 ) does not approach 
E^ when g approaches zero means that the energy when the particles are infinitely appart, 
E in f, has to be corrected. As we commented before, the integration in wavenumber diverges 
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for both (E 2 ) and for E% n f. The correlation energy is finite because both (E 2 ) and Ei n f 
diverge in the same way. This means that we expect that A will be independent of the 
choice of k max (for k ^> kp). In fact, varying the cut-off around k max we show that the 
correlation energy for h 7^ has a strong dependence on k max . For the same range of cut-off, 
A (h — > 0) is constant [Tj|. This shows again that the analysis of quantum effects, even 
in a classical plasma, is very important but that they also need to be included also in the 
energy of the particles at infinity. 



IV. ELECTROMAGNETIC FLUCTUATIONS EFFECTS IN REACTION RATES 

It is known that nuclear reaction rates in plasma proceed at a higher rate than inferred 
from laboratory experiments. This is because in a plasma, ions are not naked, but are 
surrounded by electrons that form a shielding cloud around them. The polarization clouds 
partially screen their charges resulting into a lower Coulomb barrier between them and 



an enhanced tunneling probability. The first to discuss this effect was Salpeter |jI9|| . He 
assumed adiabaticity during the collisions, which means that the electrons and ions rearrange 
themselves so the polarization cloud follows the colliding particles at all times. The reaction 
rate, R, is then found to be increased by the factor F = exp(Uo/T) |[20|| . (In the static 
case in the weak screening limit, the enhancing factor is exp(ZiZ 2 e 2 j (T Rd)), the Salpeter 
factor.) 

The exact effect of the plasma on screening has, in the last years, been a source of 
a great controversy. The solar neutrino problem has lead to a careful study of plasma 
effects in reaction rates. The solar neutrino problem is the evidence that measurements of 
electrons neutrinos fluxes in four different experiments (SAGE (Russian-American gallium 
experiment) pifl , Kamiokande p^fl , Gallex [F23] and Homestake ||24j|), showed that there is 
a deficit of neutrinos measured compared to the number predicted by the standard solar 
model. There have been a number of works p5|-pl[| discussing if the screening is in fact 
described by the Salpeter factor adressing the possibility of dynamic corrections. Carraro 
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et al. [p5fl , for example, argued that the enhancement factor should account for the velocity 



of the colliding particles. Gruzinov |29j argued that there are no dynamic corrections in 



statistical equilibrium. More recently, Bahcall et al. []32| revisited these works pointing out 
that the screening is given by the Salpeter factor. The Salpeter result is predicted by the 
Debye-Hiickel theory (with fi — > 0). Here we comment on the work of Opher and Opher 
(2000b) |3]J that examined the role of highly damped modes and how the enhancement ratio 
is modified when the limit h — > is not taken, in a classical plasma. The electromagnetic 
fluctuations in the plasma act as a background to the interacting ions. They used the 
same description as in the last Section, a mixture of a classical description of the dielectric 
permittivity, but not assuming that % — > 0, therefore quantizing the fields. The enhancement 
factor can be written as a function of the background fluctuating potential, 

F = l + Z 1 Z 2 e 2 (0 2 )[3 2 + ... = exp(z 1 Z 2 e 2 (0 2 )f3 2 ) , (8) 

where (3 — 1/T {k B = 1). (cj) 2 ) is related to the fluctuations of the longitudinal electric 
field in a plasma: (<f) 2 ) k = (E 2 ) k /k 2 . When ((f) 2 ) ^ k is integrated for all wavenumbers k, 
but only for frequencies up to ~ u pe (<^C T), the Salpeter potential is recovered. Through 
the fluctuation-dissipation theorem they evaluated the screening potential. This guarantees 
the inclusion of the normal and the highly damped modes. Because the reaction rates are 
exponentially dependent on the screening potential, any tiny changes in the potential will 
affect drastically reactions with large values of ZiZ 2 such as p — 7 Be and p — 14 N. In 
Fig. 4, the potential (4> 2 ) w fc/(0 2 )fc w= o vs UJ / UJ pe f° r the reaction p — 14 N, with T = 10 7 K 
and n = 10 24 cm -3 is plotted showing the importance of the highly damped modes. The 
solid curve is for k = 100 ko, the dashed curve for k = 10 kp and the dotted curve for 
k = kn- We see that for k = kn, the potential ((p 2 )uj,k has non-negligible values only at 
low frequencies. However, for k = 10 k D and k = 100 kjj, the potential (pk^ spreads out to 
high frequencies. For these high values of k, instead of dropping abruptly at ~ u pe , the 
curves decrease slowly with frequency. There is thus a substantial contribution for u > u pe . 
Therefore, for large wavenumbers, (i.e., k > ko) fluctuations make a substantial contribution 
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at frequencies higher than the plasma frequency. Opher and Opher (2000b) |H) find that for 
the plasmas studied, electron-proton plasmas with temperatures 10 7 — 10 s K and densities 
10 23 — 10 28 cm -3 , the reaction rates p — p, 3 He — 3 He, and p — 1A N are increased up to 
8 — 13%. Even if the exact role of quantum correction has to be more carefully examined, 
this work shows the importance of quantum corrections even for a classical plasma. 

Besides this work, another interesting aspect of the fluctuations that could affect the 
screening is short time fluctuations. In the next Section we present the result from the 
sheet model (valid in the limit of % — > 0), that for long time average reproduced the Debye- 
Hiickel prediction Eq. (7), that lead to the Salpeter screening. However, for short time 
periods, the potential reaches values much higher than the Debye-Hiickel average. In fact, 
this is the same claim made by Shaviv and Shaviv |27j that only the long time average 
of the electrostatic potential energy of a test particle in thermal equilibrium is given by 
Salpeter. However, because the reaction rates are instantaneous phenomena, the ions feel 
the oscillating potential rather than the long time average potential. All those attempts 
point out the great importance of a better understanding of the role of thermal fluctuations 
and highly damped modes. 

V. PRESENT WORK: NUMERICAL ANALYSIS WITH THE SHEET MODEL 
AND MULTIPARTICLE QUANTUM CODE 

Presently with the multiparticle quantum code and with the sheet model we are investigat- 
ing the fluctuations in the longitudinal electric field. The electrostatic sheet model is one of 
the oldest and most versatile one-dimensional models for the investigation of plasmas ]1,14 



It was used first to investigate electrostatic effects. Later, it was modified to allow motions 



on the plane of the sheets, as well as perpendicular to them ||14|| . The one species model 
|1|] considers a plasma composed of a large number of identical charge sheets embedded in 
a fixed uniform neutralizing background. The sheets are constrained to be perpendicular 
to, for example, the x-axis, and are allowed to move freely in the x direction. They also 
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are allowed to pass freely through each other. The two species model |14| consists of both 
positively and negatively charged sheets. 

Dean E. Dauger (in collaboration with Viktor K. Decyk and J. M. Dawson) developed 
a multiparticle quantum particle-in-cell code. This code models many particle quantum 
systems by combining the semiclassical approximation of Feynman path integrals with a 
parallel plasma particle-in-cell code. Work done so far on the quantum particle-in-cell shows 
that the primary computational cost of the method is in tracing the classical trajectories 
through electromagnetic fields. In the implementation, this is handled by a modified particle 
pusher and particle manager derived from a plasma particle-in-cell code. Combining these 
with a wavefunction depositor and a plasma codes field solver forms the primary components 
of the quantum particle-in-cell code. The wavefunction at t+ At is calculated using the action 
from the classical path and depositing the virtual classical particle contribution on a grid in 
real space. 

Fig. 6 show the preliminary results with the multiparticle quantum code, for the longi- 
tudinal spectra {E L ) u /{E 2 L ) w=a normalized vs. oj/u pe with 128 quantum particles (summed 
up to 10 6 classic trajectories) and with kjj = 250kBrogiie- The full curve is for k = 0.3 k D , 
the dotted curve, for k = 39 k D , the dashed for A: = 59 k D , the long dashed for k = 79 ko 
and the dash-dot for k = 99 kr>. This result agrees with the behavior of the theoretical 
prediction, displayed in Fig. 1, which shows as k increases, an increase of the intensity of 
the longitudinal spectra and the presence of the Landau damping. The presence of the non- 
normal modes is shown at high frequencies for uo > uo pe . The sheet model presents identical 
results. Fig. 6 shows for the sheet model, the longitudinal spectra {Ej) k /nk B T vs k/kr>- 
The Debye-Hiickel theory predicts the result given by Eq. (5) that for an electron plasma 
with ei(k,0) = (k 2 + k 2 D )/k 2 , gives {E 2 L ) k /nk B T = l/((k/k D ) 2 + 1). Both the temporal 
average and the theoretical prediction are plotted as solid lines, being identical. It can be 
seen that even though the long time average is the Debye-Hiickel theory, the instant value 
of the longitudinal energy oscillates and is very different from the Debye-Hiickel value. This 
points out that the short-time screening of nuclear reaction rates can be very different from 
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Salpeter (given by the Debye-Hiickel theory). 



VI. CONCLUSIONS 

We discussed in this article the effect of fluctuations in energy and in reaction rates 
showing the great contribution of highly damped modes even for the standard classical 
calculation, the Debye-Hiickel theory. Recent works showed that the highly damped modes 
are responsible for an extra energy in the tranverse component, at low frequencies. It 
also indicates the importance of quantum effects even for classical plasmas. The spectra of 
fluctuations is significantly reduced at frequencies from 0.5 — 2.0u pea k. It also pointed out the 
effect of fluctuations and highly damped modes in reaction rates. The results obtained from 
the numerical simulations with quantum and classical codes, were presented showing the 
presence of highly damped modes. The sheet model reproduced on average the Debye-Hiickel 
prediction for the level of longitudinal fluctuations. However, it shows that large deviations 
can occur for short periods of time. This indicate that the Salpeter enhancement factor can 
be correct only for long time average. All these studies show the rich and varied phenomena 
that can be studied with fluctuations in plasma. The study with the multiparticle quantum 
code will be extended to analyse carefully quantum effects in reaction rates. This code is a 
cutting edge one since there is no particle code actually existing to simulate quantum effects 
in plasma. Moreover, all these effects can have major consequences in stellar evolution 
and in the early universe. In the early universe, a modifications of the equation of state and 
reaction rates can have several consequences. It will modify the expansion rate and affect the 
production of light elements at the epoch of Primordial Nucleosynthesis, for example. These 
modifications will also affect stellar evolution. Stellar models are constructed by solving the 
basic stellar structure equations. The solution of stellar structure equations requires that 
the opacity, nuclear reaction rates, and equation of state be specified. Even tiny changes 
will affect drastically the stellar evolution and the determination of the age of the universe 
though globular clusters, for example |53 . 
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FIGURES 



FIG. 1. 

The longitudinal spectra (El) w / {E"j) u=Q vs. uo/uo pe for diferent k's. The parameters are 
T = 10 5 K, and n = 10 18 cm -3 and peak of the blackbody is 232u pe . The full curve is for 
k = O.Olkn , the dotted for 0.1k D , the dashed for 10k D and the long-dash is for k = 100k D . 

FIG. 2. 

The spectrum {Ef) ^ / {E\) W=Q vs. oj/uj pea k f° r k = k max , in the cases when with h ^ 
(full curve) and with % ^ (dashed curve). 

FIG. 3. 

The longitudinal energy, (Ej)^ , as a function of the maximum frequency of the inte- 
gration, (already integrated in wavenumber). The full curve is for % — > and the dashed 
curve is for % ^ 0. 

FIG. 4. 

The potential (0 2 )^ k / (4> 2 )k w=o vs ^/^pe fc> r the reaction p — U N, with T = 10 7 K and 
n = 10 24 cm -3 . The solid curve is for k = 100 k D) the dashed curve for k — 10 ko and the 
dotted curve for k = kp. 

FIG. 5. 

The longitudinal spectra (Ef)^/ (E^) Q vs. uo/uo pe . The full curve is for k = 0.3 ko, the 
dotted curve, for k = 39 hp, the dashed for k = 59 &d, the long dashed for k = 79 kp and 
the dash-dot for A; = 99 k D . 

FIG. 6. 

The longitudinal spectra (E\) k /nk B T vs fc/Zco 
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